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We investigate the electron momentum distribution function (EMD) in a weakly doped two- 
dimensional quantum antiferromagnet (AFM) as described by the t-J model. Our analytical results 
for a single hole in an AFM based on the self-consistent Born approximation (SCBA) indicate an 
anomalous momentum dependence of EMD showing 'hole pockets' coexisting with a signature of 
an emerging large Fermi surface. The position of the incipient Fermi surface and the structure 
of the EMD is determined by the momentum of the ground state. Our analysis shows that this 
result remains robust in the presence of next-nearest neighbor hopping terms in the model. Exact 
diagonalization results for small clusters are with the SCBA reproduced quantitatively. 
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One of the most intriguing questions concerning the 
superconducting cuprates is the existence and the char- 
acter of the Fermi surface (FS), in particular in their 
underdoped regime. This problem has been intensively 
studied experimentally with the angle-resolved photoe- 
mission spectroscopy (ARPES) There have been 

also several theoretical investigations of this problem, us- 
ing the exact diagonalization (ED) of small clusters [^-0, 
string calculations [|j , slave-boson theory ||] and the high 
temperature expansion 0]. While a consensus has been 
reached about the existence of a large Fermi surface in 
the optimum-doped and overdoped materials, in the in- 
terpretation of ARPES experiments for the underdoped 
cuprates |3j the issue of the debate is (i) why are experi- 
ments more consistent with the existence of parts of large 
FS, i.e., rather Fermi arcs or Fermi patches [p|pT| than 
with a 'hole pocket' type small FS, predicted by several 
theoretical methods based on the existence of AFM long 
range order in cuprates, (ii) how does a partial FS even- 
tually evolve with doping into a large closed one. 

The electron momentum distribution function nk = 
(*kol Z^cr 4,<TCk,<T|*ko) is the key quantity for resolving 
the problem of the Fermi surface. In this paper we study 
the EMD for |5'ko) which represents a weakly doped 
AFM, i.e, it is the ground state (GS) wave function of 
a planar AFM with one hole and the GS wave vector ko. 
In the present work we investigate the low-energy physics 
of the Cu02 planes in cuprates within the framework of 
the standard t-J model with nearest-neighbor hopping 
tii' = t and the AFM exchange J. In order to come 
closer to the realistic situation in cuprates the model is 
extended with the next-nearest-neighbor hopping tai = t' 
and the third-neighbor hopping terms tui = t" , for ii' 
representing next-nearest-neighbors and third-neighbors, 
respectively , 
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a {'^i,<y) electron creation (annihilation) operators 
acting in a space forbidding double occupancy on the 
same site. The effect of double occupancy on nk is 
not studied in the present framework of the t-J model. 
Sf are spin operators. For convenience we treat the 
anisotropy 7 as a free parameter, with 7 = in the 
Ising limit, and 7 — s- 1 in the Heisenberg model. Re- 
cent studies of the t-J model with t' , t" terms included 
have shown a very good agreement of the calculated 
quasiparticle (QP) dispersion with experimental results 
of ARPES ||l^ whereby quantitative differences between 
different Cu compounds have been attributed to different 
values of i' and t" 

Our analytical approach is based on a spinless fermion 
- Schwinger boson representation of the t-J Hamiltonian 
]l5[ and on the SCBA for calculating both the Green's 
function [p"5|-p^ and the corresponding wave function 
|l8| , ^ . The method is known to be successful in deter- 
mining spectral and other properties of the QP. In con- 
trast to other methods the SCBA is expected to describe 
the long-wavelength physics since it is determined by the 
linear dispersion of spin waves. The short- wavelength 
properties can be studied with various methods, here we 
compare the SCBA results with the corresponding ED, 
as shown further-on. 

In the SCBA fermion operators are decoupled into 
hole and pseudo spin - local boson operators: c,;.t — h] 



Ci,x = h\Sf ^ hjoi and 6^4 = h], Ci_, 
i belonging to A- and i?-sublattice, 
effective Hamiltonian emerges 

H = J2eihlhu + J2u;^. 



respectively. 
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where is the creation operator for a (spinless) hole in 
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a Bloch state with a dispersion = 4t' cos kx cos ky + 
2t" {cos2kx —cos2ky). The AFM boson operator ajj cre- 
ates an AFM magnon with the energy cjq, Mkq is the 
fermion-magnon couphng and N is the number of lattice 
sites. 

We calculate the Green's function for a hole Gk(ij-') 
within the SCBA |15-17|. This approximation amounts 
to the summation of non-crossing diagrams to all orders 
and the corresponding ground state wave function with 
momentum ko and energy e^^ is represented as 

qi 

... + ^-n/2 ^ A/kq,Gk,(u;i)...Mk„_,q„ X 

... " |0). (3) 



qi,--,qn 



Here km = kg - qi - ... -q™, iOm = Cko - Wq^ - ... -tjq,„ 
and is the QP spectral weight. The wave function is 
properly normalized (^kol^ko) — 1 provided the number 
of magnon terms n — > oo JT9] . 

The wave function Eq. (^corresponds to the projected 
space of the model Eq. (m and therefore the EMD is 
rik = (^'kol"k|*ko) = (*kol%l*ko) with the projected 
electron number operator rik = X^o- "^k (jCk.o-- Consistent 
with the SCBA approach, we decouple the latter into 
hole and magnon operators, 



rik 



X (vtj [1 + (1 - % )] + Vij (4 + a, ) j , (4) 



where r]f^ = (1 ± e~<^ (^'"^^))/2 with Q = (7r,7r). Local 

a| are further expressed with proper magnon operators 
a^. It should be noted that rik should obey the sum rule 



nk = 1 — Ch and the constraint rik < 1 
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Ch, 



where Ch is the concentration of holes . 

In general the expectation value rik for a single hole 
has to be calculated numerically and has the following 
structure 



"■k — 1 - ■r.^ko('5kko + ^kko+q) 



1 

N 
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Here the second term proportional to (5-functions corre- 
sponds to 'hole pockets'. Note that (5nk, for the case of 
a single hole fulfills the sum rule -^X^k^'^k = ^ko — 1 
and Suk < 1 . The introduction of Sni^ is convenient as it 
allows the comparison of results obtained with different 
methods and on clusters of different size N. 

For the case of Ising limit, 7 = 0, the Green's func- 
tion in the SCBA is independent of k, Gk{uj) — Go(cl'). 
Therefore it is possible to express all required matrix ele- 
ments of nk analytically and to perform a summation of 



corresponding non-crossing contributions to any order n, 
similar to Ref. The result for J/t = 0.3 (as relevant 
to cuprates) is for some selected directions in the BZ pre- 
sented in Fig. 1. We have also checked the convergence of 
iJrik with the number of magnon lines, n. For J/t > 0.3 
we find for all k that the contribution of terms n > 3 
amounts to less than few percent. This is in agreement 
with the convergence of the norm of the wave function, 
which is even faster In Fig. 2(a) we present this 5n\^ 
for the whole BZ. Here we note one interesting feature 
in the Ising limit, i.e., the dip of i5nk in the center of the 
BZ at /c ^ in agreement with Ref. 
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FIG. 1. (Srik along selected paths in the BZ for 7 = 0, 

0. 9, 0.999. Full diamonds represent ED results A^(nk — 1) for 

= 32 of Chernyshev et al. Q. Open circles represent the 
SCBA result for Srn,„ - \NZy,„. 

Now we turn to the Heisenberg model, 7^1. Here the 
important ingredient is the gap-less magnons with linear 
dispersion and a more complex ground state of the planar 
AFM. Gk(w) and Ck become strongly k-dependent. As a 
consequence nk is now in general dependent both on k 
and kg. The ground state is for the t-J model fourfold 
degenerate and we choose ko — (f, f ). Results should be 
averaged over all four possible ground state momenta if 
compared with, e.g., high temperature expansion results 
]lOt or discussed in connection with ARPES data. 

Let us first discuss the result in the limit t/J 0, 

1. e., for a static hole. In the linearized model, Eq. (j^), 
-^'^kq — 0. Therefore the hole is not coupled to the AFM 
(Zko = 1) and (5nk should be zero. However, a straight- 
forward use of Eq. (Q) leads also to non-vanishing and 
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momentum dependent Jrik = no(k). We attribute this 
momentum dependence to the improper decomposition 
of Ci^CT into Unearized pseudo spin operators instead of 
Schwinger bosons obeying the local constraints In 
the results for Sni^ presented in this paper (also at finite 
J/t) the contribution no(k) is not included. 

In Fig. 1 we present Jrik for the t-J model with (al- 
most) isotropic Heisenberg exchange, 7 = 0.999. Numer- 
ical calculations are henceforth performed for J/t — 0.3 
and N corresponding to a 64 x 64 sites cluster In 
evaluating the matrix elements we take into account only 
terms with up to n = 3 magnon lines in l^'ko) The 
sum rule of our numerical results for Sn^ is nevertheless 
fulfilled > 96%. 

Also included in Fig. 1 are results obtained with ED of 
the = 32 sites cluster Q . These data are scaled as the 
quantity N{nk — 1) which should be directly compared 
with the SCBA result Sn-k^. At momenta k = ko, ko + Q, 
however, one has to take into account contributions from 
'hole pocket' terms proportional to (5kko Eind with the 
scaling c>c A^. Thus ED data should be compared at these 
points with i5nko — \NZ\^g calculated from the SCBA. 
Note also, that the SCBA result for 7 = 0.9 (also pre- 
sented in Fig. 1) represents an intermediate step between 
the Ising limit and the Heisenberg limit: the dip at the 
r point, which is in the Ising limit well pronounced here 
disappears but the difference for directions k || ko and 
k _L ko is not yet developed. In Fig. 2(b) we present 5n\^ 
for the Heisenberg limit (7 — 0.99) in the entire BZ. In 
comparison with the Ising limit, Fig. 2(a), 5n\^ exhibits 
a very strong momentum dependence around ikg. 

The comparison of the SCBA with ED results shows 
a quantitative agreement at all points in the BZ. How- 
ever, the SCBA result is symmetric around F point in 
the direction k || ko, while small system results show a 
weak asymmetry for k = ±ko, respectively. From our 
analysis of the SCBA results for N ^ 00 and long range 
AFM spin background it follows that rik is in the ther- 
modynamic limit Ch ^ symmetric. The asymmetry 
is in Ref. attributed to the opening of the gap in the 
magnon spectrum at q ~ Q in finite systems. Within the 
SCBA the asymmetry also appears if the EMD is evalu- 
ated with ko displaced from (-1, -1) by a small amount (5ko 
(not shown here). A common feature of finite clusters is 
a non- vanishing expectation value of the current operator 
for the allowed GS wave vector. The GS with vanishing 
current may be reached by the method of twisted bound- 
ary conditions p^ , resulting in the GS momentum dis- 
placed away from The asymmetry of rik found in 
small clusters can thus be attributed to this displacement 
and is a finite size effect. In the thermodynamic limit in 
the system with AFM order the GS momentum would 
coincide with (f , f ) and no asymmetry is expected in 
nk. 

To get more insight into the structure of (5nk, we sim- 
plify the wave function, Eq. (ph , by keeping only the one- 



magnon contributions. The leading order contribution to 
(5nk is 

(5n|^^^=-ZkoMkoqGko(eko-Wq)[2'Uq-|~MkoqGko(eko-^q)] 

~-8Z^„J^(l + Zk„^), q^O, (6) 

^q 

with q = k— ko (or k— ko — Q) and v — t(sm k^x, sin fcoy). 
The momentum dependence of EMD, contained in 5rL"^\ 
essentially captures well the full numerical solution for 
the isotropic case, Fig. 2(b), as well as in the Ising limit, 
Fig. 2(a). 



\^ (a) 




FIG. 2. 5nk for (a) the Ising model 7 = 0, (b) 7 = 0.99 and 
t'^t"^0, and (c) t-t'-J model with t' = 0.25t and 7 = 0.99. 

A surprising observation is that the EMD exhibits in 
the extreme Heisenberg limit for momenta k ^ ko, ko + Q 
a discontinuity ^ Z]s_gN^^^ and Sn^"^ cx — (1+ sign qx)/qx- 
These discontinuities are clearly seen in Fig. 1, Fig. 2(b) 
and are consistent with ED results. Fig. 1. One can in- 
terpret this result as an indication of an emerging large 
Fermi surface at k ^ ikg. The discontinuity appears 
only as points ±ko, not lines in the BZ. Note, however, 
that this result is obtained in the extreme low doping 
limit, i.e., Ch — l/N and it is not straightforward to gen- 
eralize it to the finite doping regime. In the limit 7—^1 
this term does not strictly obey the constraint Sni^ < 1, 
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although due to the symmetry it does not violate the 
EMD sum rule. The singularity is weak and on introduc- 
ing a slight anisotropy, e.g. 7 < 0.999, the constraint is 
not violated. 

In Fig. 3 we present the results for the t-t'-t"-J model. 
First we introduce positive next-nearest neighbor hopping 
matrix elements t' = i/4, t" = 0, claimed to be appro- 
priate to electron doped systems such as Nd2Ce cuprates 
p3[ . The GS is now twofold degenerate, with the mo- 
menta at corners of the AFM zone, e.g., ko — (7r,0), 
with an enhanced pole residue Z\^^ = 0.54. The result 
is presented in Fig. 2(c) for the entire BZ. The discon- 
tinuity in this case disappears due to the symmetry as 
evident from v = in the leading order approximation, 
Eq. d). The effect of negative t' = -i/4, t" = is rela- 
tively weak: the GS momentum remains at kg — (-j, ^), 
(5nk at k = Q is lower than the t' = t" = result while 
the discontinuity is smaller, because Zk(, — 0.25 here. In 
Fig. 3 we additionally present results for N{ni^ — 1) ob- 
tained from exact diagonalization of a iV = V20 x \/20 
cluster. The t' — t" — results are in agreement with 
those of Ref. All ED results quantitatively confirm 
the SCBA values. A possible set of parameters appro- 
priate for reproducing the dispersion from experimental 
ARPES data is t' = -t/A, t" = t/b @. Our SCBA 
result presented in Fig. 3 is qualitatively similar to other 
t' < results. The main difference is a more pronounced 
step at k = ±ko. 
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FIG. 3. 5nk for different t' and t" . Symbols and lines repre- 
sent ED for A'' = 20 sites and SCBA results, respectively. Full 
circle and full line denote t' = t" = 0, full diamonds and long 
dashed line t' = t/4, t" = and empty diamonds, dot-dashed 
line t' = ~t/4, t" = 0. Short dashed line represents the SCBA 
result for t'^~t/4, t"^t/5. In all cases 7 = 0.999. 

In the present work we considered the electron mo- 
mentum distribution function in underdoped cuprates. 
The results of the two methods, the self consistent Born 
approximation and the exact diagonalization agree quan- 
titatively. Our analysis shows that the presence of next- 
nearest neighbor terms changes EMD only quantitatively 
if the ground state is at (f , f ) and qualitatively for suffi- 
ciently large t" > where the GS momentum is at (tt, 0). 

The main observation is however the coexistence of 



two apparently contradicting Fermi-surface scenarios in 
EMD of a single hole in an AFM. (i) On one hand, the i5- 
function contributions in Eq.(5) seem to indicate that at 
finite doping a delta- function might develop into small 
Fermi surface, i.e., a hole pocket, provided that AFM 
long range order persists, (ii) A novel feature is that also 
Suk is singular in a particular way, i.e., it shows a discon- 
tinuity at k = ko with a strong asymmetry with respect 
to kg. It is therefore more consistent with infinitesimally 
short arc (point) of an emerging large FS. For finite dop- 
ing the discontinuity could possibly extend into such a 
finite arc (not closed) FS. Note that as long-range AFM 
order is destroyed by doping, 'hole-pocket' contributions 
should disappear while the singularity in Sni^ could per- 
sist. 

Making contact with ARPES experiments we should 
note that ARPES measures the imaginary part of the 
electron Green's function. We must note that using these 
experiments in underdoped cuprates nt can be only qual- 
itatively discussed since the latter is extracted only from 
rather restricted frequency window below the chemical 
potential. Nevertheless our results are not consistent 
with a small hole-pocket FS (at least only a part of pre- 
sumable closed FS is visible), but rather with partially 
developed arcs resulting in FS which is just a set of dis- 
connected segments at low temperature collapsing to the 
point Q. The SCBA results for singular Sn]^ seem to 
allow for such a scenario. It should also be stressed that 
the SCBA approach is based on the AFM long-range or- 
der, still we do not expect that finite but longer-range 
AFM correlations would entirely change our conclusions. 
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